INVARIANT ALGEBRAIC SLICING OF THE SPACETIME* 



C. BONA and J. STELA 
Departament de Fisica, Universitat de les Illes Balears 
E-01011 Palma de Mallorca, SPAIN 

and 

J. MASSO and E. SEIDEL 
National Center for Supercomputer Applications, 605 East Springfield Avenue 

Champaign, IL 61280, USA 



ABSTRACT 

ON Using the momentum constraint, the standard evolution system is written in a 

fully first order form. The class of first order invariant algebraic slicing conditions 
is considered. The full set of characteristic fields is explicitly given. Characteristic 
^ speeds associated to the gauge dependent eigenfields (gauge speeds) are related to 

light speed. 

o 

1. Introduction 

It is well known that Einstein's field equations can be decomposed into two sets: 
O the evolution system and the constraints. There is no ambiguity in defining what 

we call constraint equations: this means equations in which there are no second 
■^j- time derivatives. Energy and momentum constraints are easily singled out in that 

way. 

The term evolution system refers to the set of remaining equations. This is 
^ just a matter of choice, because an evolution equation plus a constraint leads to 

£^ another evolution equation with the same physical solutions (the ones obtained from 

initial data which satisfy the constraints). The standard choice is to take the space 
components of the Ricci tensor, but one could choose instead the space components 
of the Einstein tensor or any other combination obtained by using the constraints 
in a suitable way, as we will actually do in the present work. 

As it is well known, the freedom to choose the cordinate gauge allows one to 
complete the evolution system in many different ways and this can lead to many 
different systems of equations, each one with its own structure. Nevertheless, we 
know that the physical solutions belong to the solution space intersection of all these 
systems so that they are physically equivalent. However, this is not of much help 
in Numerical Relativity, where hyperbolic systems play an outstanding role. They 
allow one to apply the powerful methods of Computational Fluid Dynamics: there 
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is no need for building (and testing) new 'ad hoc' codes, one can just apply what it 
is known to work and concentrate himself in the physics of the problem considered. 

In this work we will look for hyperbolic evolution systems for a wide class of co- 
ordinate gauges. We shall consider first order slicing conditions which are invariant 
under (time dependent) transformations of the space coordinates. We shall take for 
simplicity a vacuum spacetime and we will use a zero shift vector, but most of the 
results apply also to the general case. 

2. Evolution Equations 

The standard 3+1 evolution system is given in the vacuum case by: 

dt7,-j = -2a Kij (1) 

dtKij = -a,- ;j + a [4 3) " 2tf J + tr K K tl ] (2) 

This well-known system is of first order in time, but second order in space. To 
obtain a system which is also of first order in space, we will introduce auxiliary 
variables which correspond to the space derivatives, 

A k = d k lna } D klJ = 1/2 d k ^ l3 . (3) 

One could then simply insert these quantites into the standard ADM equations to 
obtain a first order system. However, doing so blindly does not have any particular 
advantage. As we show, a careful choice of variables will transform the equations 
into a special flux conservative, hyperbolic form that is especially suited to numerical 
treatment. 

In particular, the evolution system (2) can then be written as 

3 t K t] + d r (X^) = a S t] . (4) 
where the terms Af- are given by 

A£ = D% + 1/2 b\ (Aj + 2V,- D; r ) + 1/2 8) (A, + 2V t - DJ), (5) 
and we have noted for short 

V k = D k [ - D r rk . (6) 
Sij is a source term involving only the fields themselves and not their derivatives: 

- (2D kr k - A r )(D tjr + D jir ) + A t (Vj - 1/2 D ]k ) + A 3 (V t - 1/2 D tk k ) . (7) 

It is clear that one needs also to evolve the space derivatives. The simplest way 
of doing so is just to take the time derivative of (3) and interchange the order of 
space and time derivatives: 

d t A k + d k (aQ) = 
d t D kij + d k (a K^) = (8) 



where we have noted 
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a Q , 



(9) 



so that the choice of Q will determine the slicing. 

The vector Vk is a very interesting quantity. One can compute its time derivative 
from (8) but we will use the momentum constraint to transform this equation into 



8 t V k = aA r (K\ - tr K 6 r k ) + a [(D k s r - 2D r s k ) - 6 s k (D r j - 2D^ r )} K\ . (10) 



We will consider Vk as independent quantities to be evolved with equation (10). In 
that way, (6) becomes an algebraic constraint between V and the spatial metric 
derivatives D (the algebraic form of the momentum constraint). As we will see, this 
is crucial to ensure the hyperbolicity of the evolution system. 

This set of equations (1,4,8-10) has the special form we seek. The entire non- 
linear system of the Einstein evolution equations is now written in the form of first 
order balance laws as 



which is familiar from many branches of physics. The structure of this system can 
be now easily investigated. 

3. Invariant algebraic slicing 

Note that we have not yet specified the slicing, because the proper time deriva- 
tive Q of the lapse function is yet to be given. We are interested in invariant slicing 
conditions. This means that the spacetime slicing provided by our coordinate con- 
dition must be invariant under any transformation of the space coordinates of every 
slice. We must use then slicing scalars, like a, Q or tr K and their proper time 
derivatives (note that the shift "vector" does not behave as a slicing vector; it is a 
vector under time independent transformations only). 

We also want to use an algebraic condition. If we restrict ourselves to zero order 
scalars, we can play only with a and we get either a geodesic slicing or one of its 
generalizations. If we allow also first order scalars, we get both Q and tr K into the 
play. The most general homogeneous algebraic condition is then 



where / is an arbitrary function. 

The geodesic slicing is then included as a subcase with / = 0. The maximal 
slicing is included also as a limiting case when / diverges. The / = 1 case cor- 
responds to the harmonic slicing. Another interesting case is the '1+log' slicing, 
obtained when / = 1/ct; it mimics maximal slicing near a singularity, when the 
lapse collapses to zero. The term '1+log' arises from the expression of a in terms of 

one obtains when integrating (12) in the eulerian case (zero shift vector). Note 
however that the invariance of (12) ensures that one can apply it to obtain the same 
slicing even with a nonzero shift vector. 







4. The characteristic fields 



Although the system of evolution equations (1,4, 8-10) is first order and flux 
conservative, we have not examined its hyperbolicity. When dealing with a first 
order system, one does not discuss hyperbolicity in general. One must first choose 
a fixed space direction (we will take the the k coordinate axis) to discuss hyperbol- 
icity by considering only space derivatives along the selected direction. This is a 
framework which does not match the usual one for second order equations, where 
there is no need for choosing a priori a direction and all derivatives are dealt with 
at a time. 

It is clear that the normal lines will be characteristic (with zero speed). The 
corresponding characteristic fields are the lapse and the metric components, the 
combination appearing in (12), the vector V and 

A k -fD k r r , A k ,, D k 'ij (k'jtk). (13) 

Also, the light cones are characteristic surfaces (speeds ictyj r y kk ). The correspond- 
ing eigenfields are 

\fy^ K lk i =p X k lk ,, (14) 
One gets also characteristic cones associated to the eigenfields 

yf^trKTy/f X (15) 

with characteristic speeds 

±<V/ 7 fcfc (16) 

and we will call these 'gauge speeds' because their explicit dependence on the slicing 
condition. Gauge speed coincides with light speed only in the harmonic case. It 
becomes infinite for a maximal slicing, which can be considered as a limiting case 
of our condition (12). 

It is clear that negative values of / will lead to imaginary gauge speeds. More- 
over, the set of eigenfields is complete only if / ^ 0. This means that the evolution 
system will be strictly hyperbolic iff / > 0. Note also that gauges with f < 1 
will have poor singularity avoiding behaviour because gauge speed would be lower 
than light speed. Therefore, cases with / > 1 will look more appealing for most 
Numerical Relativity applications. 
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